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Abstract 

Recent works have shown that non-equihbrium concentration fluctu- 
ations arise during free difi^usion in simple, three dimensional fluids. In 
this paper, I extend the theory of non-equilibrium fluctuations to two di- 
mensional systems, and I evaluate the root mean square amplitude of the 
fluctuations in the presence of an uniform concentration gradient. I find 
that, in this case, fluctuations are expected to be extremely strong. Ex- 
perimentally, two-dimensional liquid systems are approximated by freely 
suspended liquid films; diffusion is studied with I discuss the effect of 
the drag due to the fluids surrounding the film, that is the main known 
difference between real fllms and ideal two-dimensional liquids. In or- 
der to check experimentally this theory, I provide an evaluation of the 
power spectrum of the fluctuations that should arise during the diffusion 
of a spot with a gaussian proflle, that can be directly compared with the 
results from a fluorescence recovery after photobleaching (FRAP) exper- 
iment. I propose that this kind of experiment could provide informa- 
tions on the physical properties of membranes. PACS: 68.15.-fe Liquid 
thin fllms, 68.65.-k Low-dimensional systems, 47.10.-|-g Fluid dynamics, 
05.40.-a Fluctuations phenomena, 05.70.Ln, Nonequilibrium thermody- 
namics, 66.10.Cb Diffusion. 

1 Introduction 

Pick's flow arises whenever a concentration gradient is present in a multicompo- 
nent fluid. The presence of this flux is the evidence of non-equlibrium condition; 
recently, fluctuations much stronger than equlibrium ones have been reported 
in systems made by couples of ordinary miscible fluids. [U [2j |3l |4] These non- 
equilibrium fluctuations are originated by the termally excited velocity fluctu- 
ations, that displace volumes of fluid in layers with different concentrations. 
Also, it has been shown that this effect is the only origin of the Pick's flow [S]. 

Until now, the performed experiments and the developed theory concerned 
three dimensional fluids, where fluctuations, though stronger than the equilib- 
rium ones, are always very weak, and can be considered the low-frequency tail 
of a process that happens mainly at the molecular level. In this paper, I use the 
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theory of non-equilibrium fluctuations [5] to show that, in two dimensions, the 
fluctuations are very strong, and that they have a macroscopic, non negligible 
component. The calculation of the Pick's flow in an ideal two-dimensional liq- 
uid gives the already known divergence, that here I interpret as the effect of the 
divergence of the fluctuations at long wavelengths; the usual Saffman formula 
can be obtained when considering the drag due to the surrounding fluids. 

Usually, diffusion coefficient in freely suspended liquid films is measured by 
FRAP (fluorescence recovery after photobleaching) ; in a FRAP experiment, a 
spot is bleached on a film that contains a fluorescent dye, by a strong laser radia- 
tion; then, the diffusion of the bleached fiuorophores is observed as a progressive 
smearing of the bleached spot. In this paper, I present the calculation of the 
power spectrum expected for the fluctuations arising during the diffusion of a 
component, that starts with a gaussian distribution, similar to what happens 
in FRAP experiments. The results should be valid for completely liquid films; I 
propose that the measurement of deviations from the expected spectrum could 
be used to study the effect of rigid structures in the film, or the interaction with 
surrounding media. 



2 Root mean square amplitude of the fluctua- 
tions. 

The time evolution of the concentration c {x, t) is given by convection and dif- 
fusion: 

d 

— c (x, t) = ~u (x, t) ■ \7c , t) + DV^c (f , t) (1) 

where u{x^ t) is the velocity field and D is the diffusion constant. I assume that 
the concentration gradient is a constant, since it is always much more than the 
fluctuations. In Fourier space: 

" iLoc {q, uj) ~ —it (q, u)) ■ Vcq — q^Dc {q, uj) (2) 

By solving with respect to the concentration, I obtain the correlation functions: 

, (u{q,u}) ■ zu* (q',uj') ■ z) 
c iq, c* ) . Vol \ ^ / (3) 

where z is the versor of the concentration gradient. The velocity correlation 
ftmction for a liquid is: 

{u{q,^) ■ In* (q',u') -1)^5 {q ~ q') S (uJ - Uj') ^^J^ , ^ ,1 (4) 
In both 3 and 2 dimensions [21 IS E] : 

(c (<f, 0.) C* iq', LO')) ^ 6 (q- q') 5 {u^ - u') c \ / 7 , (5) 
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It sould be noticed that, in 2d, p is the superficial mass density. 

The static power spectrum of the fluctuations can be derived, by integrating 
over Lo. The diffusion time is much longer than the viscous time, since D <^ v. 

KrT o 1 1 9^ - (9 • z) 
(c (,-, t) c i,^, t)) ^ S (,^+ <f) ^Vc^--^^ (6) 

Root mean square value of the fluctuations of concentration, in 3d, by using 
polar coordinates: 

^rms = ^^c'^l^'\^j\in{0)dej^q'dq'-^ (7) 

where Q is the wavevector of the longest wavelength fluctuations. At the begin- 
ning, the concentration of the solution increases by Ac along a length L. The 
gradient is thus Vc = Ac/L. The time needed for the fluctuations in order to 
fully develope is their diffusive time. The most intense fluctuations will be found 
approximately at the diffusive time of the macroscopic concentration proflle; at 
that time Q ~ tt/L. 
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Using Stokes-Einstein relation [9] : 



(9) 



^ - (10) 

for writing D in terms of the shear viscosity 77 and the radius a of the diffusing 
particles, I obtain: 

Crms = Ac^l^^l^ (11) 

This equation shows that fluctuations are essentially a molecular effect in 
3D: as the size L of the sample increase with respect to the molecular size a, 
the fluctuation amplitude decrease. 

Root mean square value of the fluctuations of concentration, in 2d: 

Following the previous calculation: 

2 _ KbT^, 1 



Ac^— (14) 



47r4p lyD 



3 



The term \/Q multiplied by the gradient of c gives the total concentration 
difference. 

This equation shows that, for an ideal 2D fluid, the relative fluctuation 
amplitude does not vanish as the size L of the system increases. Fluctuations are 
not a molecular effect in 2D. Moreover, it can be noticed that the fluctuations, 
as expected from Eq. , can be of the of the same order to the macroscopic 
concentration difference. We could also consider cases in which Eq. p4)) give 
unreal result Crms > Ac: this means that linearization is not allowed. 

In the case of a free-standing liquid film, the diffusion coefficient D must 
be evaluated by using Saffmann equation, that holds for a thin sheet of liquid, 
surrounded by a less viscous fluid [3 [TUl [TT] : 



where qc — 2r]i/ (hrjo), rjo is the shear vescosity of the fluid film, rji is the shear 
viscosity of the surrounding fluid, h is the thickness of the fluid sheet, a is the 
radius of the molecules, and / « f .12294 is an adimensional constant, close to 1. 
Since we are dealing with a two dimensional fluid, the mass density of the liquid 
composing it is p/h, and thus rjoh = pv. The meaning ot the critical wavelength 
qc will be discussed in more depth later. I obtain: 



For a soap film Crms ~ 0.07Ac. 

Comparing eq. pT|) and eq. (fT6|) we can notice that the greatest root mean 
square amplitude reached by the fluctuations in three dimensions decreases as 
the initial concentration distribution becomes macroscopic, with the square root 
of the size of the region where the gradient is present. In two dimensions the 
process is completely independent on the size of the region where the gradient 
is present. There is indeed a macroscopic characteristic size, that is l/qc, but 
the dependence of the fluctuation amplitude on this characteristic size is only 
logarithmic. 

3 Non-equilibrium fluctuations with a non-constant 
gradient 

In experiments on free-standing liquid films, the gradient is limited to a small 
region of the sample. For example, we can use fiuorescent molecules, and mark 
a part of the molecules by photo-bleaching, and then we can follow the diffusion 
and formation of fluctuations in fluorescence. In order to evaluate the results of 
such an experiment, we must know what happens in the case of a non-constant, 
time-dependent gradient. Here are the calculations. 




(15) 



= Ac 



1 



(16) 
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di 



c {x, t) = -u (f , t) ■ Vc (f , t) + DV^c [x, t) 



Approximation: c (x, t) = 5c {x, t) + cq {x, t) 
At zeroth-order: 

d 

—Co {x,t) = D\/^co {x,t) 

At the first order in the fluctuations: 
d 



dt 



5c (f , t) = -u {x, t) ■ Vco (x, t) + DV^(5c (f, 



Ansatz: 5c {x, t) = ^ (x, t) ■ Vcq {x, t). 



d 



(f , • Vco = -u (f , • Vco + DV^^ (x, • Vco 



(17) 



(18) 



(19) 



(20) 



There are also other mixed terms, d/dxVcQ x d/dx^, which are negligible if the 
fluctuations are on spatial scales smaller than the macroscopic concentration 
profile. Moreover, we are interested only in the component of ^ along Vc, so 
we can arbitrarily add a second equation, thus obtaining the following vector 
equation: 



d_ 
dt 



^ (x, t) = -M (x, t) + DV^^ (x, t) 



(21) 



(22) 



(23) 



In the Fourier space: 

§-/i^,t) = -uiq,t) - Dq^^{q,t) 

Solving we get: 

rn(g,t')e-^«'(*-*')dt' 

Now, we calculate the correlation matrix of ^: 

iiq,t)-i^{q',t))= J' dt'dt" {u{q,t')-u*^ (g', O) e-^l«^'(*-*')e-^l«1'(*-*") 

(24) 

The velocity correlation is a Dirac delta, whence q = q'; moreover, since u{q) is 
orthogonal to q, the matrix is the operator which projects orthogonally to q: 



i{q, t) ■i^{q',t))=5 {q- <f) {q, t) 



(25) 



S^{q,t) 



q ■ q 



S^iq,t) 



?j (g, t) = f dUt" {u {q, t') ■ u* {q', t")) e-^l«1^(2*-*'-*") 
Jo 



(26) 



(27) 
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By using the definition of spectrum in w we get: 



Su {q, 



Now we calculate the integral in dt': 



, ^ e-^^^^^ + l-2cosMe-^^^^ ^ , 



In the limit of t > l/{Dq^): 
{q, t +00) 



do; 



1 



^2 + D^q^ 



(28) 



(29) 



(30) 



Anyhow, the width of the Maxwellian curve is Dq^, which is much less than 
the width in uj of the velocity correlation; for this reason we can approximate 
as follows: 



{q, t) = 5*^ (g, cj = 0) / dcj- 
By integrating: 



-2DqH + 1 - 2cos (a;i)( 

w2 + £)2g4 



-Dq^t 



1 _ e-2^9 * 

{q,t) ^ TT ^-^ Ss{q,uj = 0) 

From Ref. [T^], we can evaluate the velocity fluctuation: 

l-e-2^?'* KbT 1 



{q,t) = TT- 



Dq2 An^rjoh q {qc + q) 



(31) 



(32) 



(33) 



Here we can also write 779/1 — pov, with the superficial density. Now we 
apply the SafFman formula for the diffusion coefficient in 2D: 



S^{q,t) 



1 1 _ e-^Dq-t 
nln(^] lHqc + q) 



(34) 



Ritorniamo a delta c, e scriviamo la trasformata del prodotto come con- 
voluzione delle trasformate: 



6c {q, t) = -i / AqcQ{q ,t)q • ^ (g - q , t) 



(35) 



Calcoliamo lo spettro; in questo caso non si fa < dc[q)dc{q') >, perche' di 
fatto c'e' una correlazione anche se q ^ q' . 



\Sc {q, on = / dq' \co {q', t)\^ q' ■ {q - q' , t) ■ q' 



(36) 
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Ora usiamo ancora la approssimazionc su c: che il suo modulo sia divcrso 
da solo per valori di g' piccoli, molto piu' piccoli dei q a cui siamo interessati. 



\5c [q, tf)= Cl {q, t) S^-{q, t) (37) 
dove 

Cl iq,t) = [ dq- |co iq-,t)f ' }f ' (38) 

rappresenta una misura della quantita' di gradiente macroscopico. Si noti che 
dipende solo dalla direzione di q, e non dal modulo. Una volta fatto lo spettro 
di potcnza dell'immaginc, che c' la somma dcgli spcttri di potcnza da co c Sc, 
possiamo valutare Cq integrando per valori di q piccoli, e poi dividere lo spettro 
risultante per Cq. 

Amplitude of the fluctuations: 



4^, (x, t) = Vci {x, t) ■ {^^ {x, t) ■ e {x, t)) ■ Vco {x, t) (39) 

From the definition of the Fourier transform: 

4ms t) = VcJ {S, t) ■ j S^{q, t) &q- Vcq (f , t) (40) 

Now we rewrite the equations in polar coordinates; we use the definition 
of 5^ and we add a low-frequency cut-off, representing the macroscopic length 
scale on which cq is present. 



c2 

rms 



2 /'27r p-\-oo 

(x,t)= Vco{x,t) J di?cos2(t?)y qdqS^{q,t) (41) 



Now we put the expression of S: 

elms ix,t) = Vco {X,t) TTT / n2 (n J- n\ ^"^^^ 

Trlnl-^) ^0 r (9C + 3) 

\aqc ) 

We change the integration variable: 



c: 



rms 



2 



{x,t)= Vco{cc,t) y ^2Dtf{ql2Dt) (43) 



TTlnf^l 

\aqc ) 

where: 

The integral can be solved in terms of special functions. 



/ [x) = -^^e-^ [-TrErfiVi + Eia;] - 7 + 2^Jt^ - logx (45) 
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where Ei is the integral exponential function, Erfi is erf{ix)/i, and 7 is the 
Euler gamma function. The function / (x) is shown in Fig. [1] Limits: for 

X ^ f {x) Ki yj - log {x) /2; for x +cx) / (x) « 1/ \fxpK 
Root mean square: 



The numeric factor (the root), is about 1/3. The term VclDt is the difference 
of concentration over a distance equal to the diffusion length. 

4 Non-equilibrium fluctuations in a FRAP ex- 
periment on a freely-suspended liquid fllm 

In the previous section, a two-dimensional fluid was analyzed, in which a con- 
stant concentration gradient was present. In this section, I analyze an actual 
system, that could be experimented. In the real world, a thin, freely suspended 
liquid layer is considered a system that follows two-dimensional hydrodynam- 
ics, on length scales longer than the layer thickness. Many experiments have 
been performed on soap films, that consist of two monolayers of amphiphilic 
molecules around a layer of water. The layer thickness ranges from the 4nm 
of the Newton black film to many micrometers; features down to a fraction of 
a millimeter should be considered two dimensional |13'. Also smectic A liquid 
crystals can form films, by a different effect. It spontaneously form a layered 
structure, with the elongated molecules alligned to the surface normal. When 
a film is drawn from this kind of liquid crystal, it is formed by a set of mono- 
layers, ranging from two to many hundreds |llj : the thickness can be selected 
with extreme accuracy. Because of the absence of translational order in the 
monolayers, the film behaves as a two dimensional system. 

Many papers [13l [Ml [15l [16l [T7| describe experiments in which turbulent 
flows have been detected, and discuss the agreemen with two dimensional hy- 
drodynamics. One result is that the role of the air surrounding the film cannot 
be neglected. In turbulent flows, it damps the turbulent vortices. From a the- 
oretical point of wiev, an exact evaluation of this effect is extremely difficult; 
for example, Rivera et al. '181 modeled this damping as a linear drag term 
in the two dimensional Navier-Stokes equations, in order to describe velocity 
fluctuations observed in a turbulent flow. They evaluated the amplitude of the 
drag term from experimental data; from this value, they found the dissipation 
due to air friction to be a significant energy dissipation mechanism in their 
experimental system. 

Another well known effect of air friction on thin liquid films can be found 
in diffusion of particles, whose diameter is about the thickness of the layer. 
Neglecting air friction, the mobility of a particle in a two dimensional system 
should be infinite, that is, a particle subjected to a steady force should accelerate 



(f , t) = Vco (x, t) Vot- 
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(46) 
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indefinitely. This would result in an infinite diffusion coefRcient. Saffman [7] 
found that air friction reduces the mobility to a finite value: the author evaluated 
the diffusion coefficient and found that it strongly depends on the viscosity of 
the surrounding fluid, and diverges as the viscosity vanishes. Saffman formula 
has been tested both in soap films [T^ and in liquid crystal films [llj . 

Diffusion has been measured in liquid crystal films, by FRAP experiments. 
The observation of the diffusion of the bleached spot sould show strong fluctu- 
ations. In the following, I describe an algorithm for evaluating the fluctuations, 
arising during the diffusion of a gaussian-shaped spot; I also include the effect 
of drag of the surrounding fluids. 

Now we evaluate Cq for a Gaussian spot. First of all, we use the circular 
symmetry: 

C2(t)= /g'dg'd^|co(q',t)|'g'2sin2(^)= / -±q'3 \c, {q' ,t)f (47) 



2 

and we put an explicit solution of the diffusion equations: 

Co(9,0 = coe-«'/'^''*' (48) 
where a'^{t) = + 2Dalt) 

1 r+oo 2 

Cl it) ^'fj^ dgg3e-.V^^(*) ^ |,4(^) (49) 
Now we evaluate the whole spectrum: 



\aqc J 

(50) 



We rescale the variables, by placing x = q/a (t) and Sc (x) = Sc (x) /cq. 

^""^ = + AnHn(^] xHqc/a{t) + x) ^^^^ 
V age / 

For 2D(Tq {t)t ^ 1, i.e. asymptotically for long times: 

2 1 1 - e"^^ 

Sc(x)Me ^ + -, ^^5- -— (52) 

AnHn ^) x^qc / <J (t) + x) ^ ^ 

The result is shown in Fig. [5] 
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5 Other notes.. 



Theoretical evaluation of diffusion coefficient in thin sheets of viscous fluid, 
surrounded by a less viscous fluid 

Diffusion in free-standing liquid films: soap films [TU], liquid crystal films 

Effect of viscosity of air in fiows; modeled as a linear drag term, with 2D 
Navier- Stokes equations [TS]. 



6 Fluctuations as the origin of Fick's flow 

and 

^ ^ lu{q,uj)- zu* [q',uj'] ■ z) 
(c (^, u* (,', ) . Vc„ \ ^ / (53) 

The power spectrum of the velocity fluctuations has been evaluated for a 
film surrounded by a fiuid [121. In particular, for = 0: 

5„ {q,u = 0) = . (54) 

A-K-^rjoh q [qc + qj 

where: 

qc = (55) 
hriQ 

Kb is the Boltzmann constant, T the temperature, etao and etai are the shear 
viscosity of the film and the surrounding fluid, and h is the thickness of the fllm. 

The power spectrum has a. q~^ dependence for q^ qc] at smaller wavevec- 
tors, the divergence saturates to g^^, due to air damping. For a soap film in 
air, the values of densities are po — lO'^Kg/m^ and pi = 1.3Kg/m'^. The ef- 
fective viscosity of the film is z/q — 1.6 • 10~^m^/s, [18], evaluated by means of 
Trapeznikov relation [19]. The viscosity of air is vi = 1.43 • 10~^m^/s [20]. From 
these values, for a 2/im thick film, qc ~ 1.2 • 10*m^^, corresponding to a 0.5mm 
wavelength. 

From the power spectrum of the velocity fluctuations, we can derive the 
diffusion coefficien[5]: 

D^TT f Su{q,OJ = 0)d^q (56) 
Ja 

where the cut off wavevector A is of the order of the inverse of a, the radius of 
the diffusing particle. By using the power spectrum of eq. (|54p . and using polar 
coordinates: 

KbT f 1 



D « / ^— dg, (57) 

^T^mh J^qc + q 

All the wave vectors less than qc contribute equally to concentration diffusion; 
shorter wave lengths, down to the molecular scale, give a decreasing contribu- 
tion. Integrating: 

D = ^\ln^-IM2)] (58) 



2'KrjQh 



In ^- In (2) 
aril 
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This is the weh known Saffman formula for the diffusion coefficient of a particle 
on a liquid film [H [lOl [U] . 

Confronto con una teoria fenomenologica [H] 



7 Introduzione vecchia 

The diffusion coefficient can be evaluated once we know the correlation func- 
tion of the thermally excited velocity fluctuations; this means that the velocity 
correlation function must be strongly affected by air friction. In this paper, I 
evaluate the correlation function of the velocity field of a thin suspended liquid 
film, by solving the hydrodynamic equation of the whole system, including the 
liquid film and the fluid surrounding it. 

Then I discuss the result, showing that both the temporal and the spatial 
statistical properties are the ones of a two dimensional, hypotetical fluid, only for 
temporal and spatial frequencies higher than critical values; at lower frequencies 
they show significative deviations due to the interactions of the liquid layer with 
the surrounding fiuid. 

Last, I derive approximately the value of diffusion coefficient, and I compare 
it with Saffman formula. 

8 Definitions 

Fourier transform: 




(59) 



Fourier anti-transform: 




(60) 



Dirac delta: 




(61) 



Convolution: 




(62) 



Convolution / product: 



F[f{x)g* m^fiq)^9i'f) 



(63) 
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